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Abstract 

In a series of papers, we investigate the reformulation of Epstein-Glaser renorma- 
lization in coordinate space, both in analytic and (Hopf) algebraic terms. This first 
article deals with analytical aspects. Some of the (historically good) reasons for the 
divorces of the Epstein-Glaser method, both from mainstream quantum field theory 
and the mathematical literature on distributions, are made plain; and overcome. 

I Introduction 

This is the first of a series of papers, the companions [[I], often being denoted respectively 

II and III. 
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We find it convenient to summarize here the aims of these papers, in reverse order. Ever 
since Kreimer perceived a Hopf algebra lurking behind the forest formula |J, the question of 
encoding the systematics of renormalization in such a structure (and the practical advantages 
therein) has been in the forefront. Connes and Kreimer were able to show, using the model 
as an example, that renormalization of quantum field theories in momentum space is encoded 
in a commutative Hopf algebra of Feynman graphs TC, and the Riemann-Hilbert problem 
with values in the group of loops on the dual of H ||, |5|. The latter makes sense only in 
the context of renormalization by dimensional regularization || fTJ, physicists' method of 
choice. Now, whereas it is plausible that the Hopf algebra approach to renormalization is 
consistent with all main renormalization methods, there is much to be learned by a systematic 
verification of this conjecture. Paper III focuses on combinatorial-geometrical aspects of this 
approach to perturbative renormalization in QFT within the framework of the Epstein- 
Glaser (EG) procedure 0. 

One can argue that all that experiments have established is (striking) agreement with 
(renormalized) momentum space integrals 0. Be that as it may, renormalization on real 
space is more intuitive, in that momentum space formulations "rather obscure the fact that 
UV divergences arise from purely short-distance phenomena" |IOfl . For the questions of 
whether and how configuration space-based methods exhibit the Hopf algebraic structure, 
the EG method was a natural candidate. It enjoys privileged rapports with external field 
theory |TI|, 12, 13], possesses a stark reputation for rigour, and does not share some limitations 
of dimensional regularization — allowing for renormalization in curved backgrounds H ] , for 
instance. 

In spite of its attractive features, EG renormalization still remains outside the mainstream 
of QFT. The (rather rigorous) QFT text by Itzykson and Zuber has only the following to 
say about it: "... the most orthodox procedure of Epstein and Glaser relies directly on the 
axioms of local field theory in configuration space. It is free of mathematically undefined 
quantities, but hides the multiplicative structure of renormalization" |15|, p. 374]. Raymond 
Stora, today the chief propagandist of the method, had commented: "In spite of its elegancy 
and accuracy this theory suffers from one defect, namely it does not yield explicit formulae 
of actual computational value" [TO. Indeed. 



Over the years, some of the awkwardness of the original formalism was dispelled in 
the work by Stora. The "splitting of distributions" was reformulated in as a typical 



problem of extension (through the boundaries of open sets) in distribution theory. Moreover, 
in [pl| it was made clear that an (easier) Euclidean analog of the EG construction does 
exist. Beyond being interesting on its own right (for instance for the renormalization group 
approach to criticality), it allows performing EG renormalization in practice by a (sort of) 
"Wick rotation" trick — the subject of paper II of this series. 

When tackling the compatibility question of EG renormalization and the Connes-Kreimer 
algebra, two main surviving difficulties are brought to light. The first is that, while the Hopf 
algebra elucidation of Bogoliubov's recursive procedure is defined graph-by-graph, in the 
EG approach it is buried under operator aspects of the time-ordered products and the S- 
matrix, not directly relevant for that question. This problem was recently addressed by 



Pinter [19, E0[ and also in [21]; the last paper, however, contains a flaw, examined in III. 



The second difficulty, uncovered in the course of the same investigation, has to do with 
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prior, analytical aspects of Epstein and Glaser's basic method of subtraction. For, it was 
curious to observe, the extension method by Epstein and Glaser has remained divorced as 
well from the literature on distributions, centered mainly on analytical continuation and 
"finite part" techniques. One scours in vain for any factual link between EG subtraction and 
the household names of mathematical distribution theory. 

And so the vision of casting all of quantum field theory in the light of distribution 
analysis |22], ^ has remained unfulfilled. 

In the present paper we are concerned with the second of the mentioned difficulties. 
This means in practice that we deal with primitively UV divergent diagrams. (Nonprimitive 
diagrams are dealt with in III.) By means of a seemingly minute departure from the letter, 
if not the spirit, of the EG original prescription we succeed to deliver its missing link to the 
standard literature on extension of distributions. Then we proceed to show the dominant 
place our improved subtraction method occupies with regard to dimensional regularization 
in configuration space; differential renormalization [24]; "natural" renormalization [25]; and 
BPHZ renormalization. 

The benefits of the improved prescription do not stop there: it goes on to remarkably 
simplify the task of constructing covariant renormalizations in II, and the Hopf-geometrical 
constructions in III. 

An important sideline of this paper is the use of the theory of Cesaro summability of 



distributions |26], \Hj in dealing with the infrared difficulties; this helps to clarify the logical 
dependence of the BPHZ procedure on the causal one, already pointed out in P8|. Improved 
BPHZ methods for massless fields ensue as well. 

The main theoretical development is found in Section 2. Afterwards, we proceed by way 
of alternating discussions and examples. In order to deliver the argument without extraneous 
complications, we work out diagrams belonging to scalar theories. Most examples are drawn 
from the massless ip\ model: masslessness is more challenging and instructive, because of the 
attendant infrared problems, and more interesting for the renormalization group calculations 
performed in III. Eventually we bring in examples in massive theories as well. 



2 Renormalization in configuration space 
2.1 The new prescription 

All derivatives in this paper, unless explicitly stated otherwise, are in the sense of distribu- 
tions. We tacitly use the translation invariance of Feynman propagators and amplitudes; in 
particular, the origin stands for the main diagonal. 

Let d denote the dimension of the coordinate space. Typically, d will be An. An unrenor- 
malized Feynman amplitude f(T), or simply /, associated to a graph T, is smooth away from 
the diagonals. We say that T is primitively divergent when f(T) is not locally integrable, but 
is integrable away from zero. Denote by < Fp r i m (R e! ) ^ L[ oc (R d \{0}) this class of amplitudes. 

By definition, a tempered distribution / e iS>'(R d ) is an extension or renormalization of / 

if 

/M :=</,$= / f(x)<P(x)d d x 
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holds whenever belongs to S(Mr \ {0}). 
Let 



f(x) = 0(\x\- a ) 



as x 



0. 



(1) 



with a an integer, and let k = a — d > 0. Then, / ^ L 1 1 oc (IR d ). But / can be regarded as a 
well-defined functional on the space iSfc + i(R d ) of Schwartz functions that vanish at the origin 
at order k + 1. Thus the simplest way to get an extension of / would appear to be standard 
Taylor series surgery: to throw away the fc-jet of <p & t the origin, in order to define / by 
transposition. Denote this jet by j^cf) and the corresponding Taylor remainder by Rq4>. We 
have by that definition 



</, <j>) = (f, Rfo] 

Using Lagrange's integral formula for the remainder: 



(2) 



R k <f>(x) 



(k + l) £ 

|/3|=fe+l 



X' 



dt(l-t) k dP<j){tx), 



where we have embraced the usual multiindex notation, and exchanging integrations, one 
appears to obtain an explicit integral formula for /: 



f(x) = (-)^ 1 (k + l 



\P\=k+l 



dt 



fk+d+ 



rf 7 



X 



(3) 



Lest the reader be worried with the precise meaning of (H), we recall that in QFT one 
usually considers a generalized homogeneity degree, the scaling degree [29]. The scaling 
degree a of a scalar distribution / at the origin of M d is defined to be 



Of = inf{ s : lim X s f(Xx) 



0} for / e S'(R d ), 



where the limit is taken in the sense of distributions. Essentially, this means that f(x) = 
0(\x\~ a f) as x — » in the Cesaro average sense |30| . Then [of] and respectively [of] — d 
— called the singular order — occupy the place of a in (P and of k. 

The trouble with ([3D is that the remainder is not a test function, so, unless the infrared 
behaviour of / is very good, we end up in @) with an undefined integral. In fact, in the 
massless theory / is an homogeneous function with an algebraic singularity, the infrared 
behaviour is pretty bad, and — d is also the critical exponent. A way to avoid the infrared 
problem is to weight the Taylor subtraction. Epstein and Glaser H introduced weight 



functions w with the properties w(0) 
4> i — * W w (p on 5(R a! ) given by 



1, w^ a \0) = for < | ct | < k, and projector maps 



[x) - w(x)j (p(x) 



(4) 



The previous ordinary Taylor surgery case corresponds to w = 1, and the identity 



W w (w<f>) = wW 1 (f> 
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tells us that W w indeed is a projector, since Ww^wx 1 ) = for |7| < k. 

Look again at @. There is a considerable amount of overkilling there. The point is 
that, in the homogeneous case, a worse singularity at the origin entails a better behaviour 
at infinity. So we can, and should, weight only the last term of the Taylor expansion. This 
leads to the definition employed in this paper, at variance with Epstein and Glaser's: 



T w <j>(x) := <j>(x) -jt\<l>K*) -*>W E ^ W (0). 

z — ' a\ 

\a\=k 



(5) 



Just w(0) = 1 is now required in principle from the weight function. 

An amazing amount of mathematical mileage stems from this simple physical observation. 
To begin with, T w is also a projector. To obtain an integral formula for it, start from 



T 

J- it 



[l-w)R k -\ 



+ wR 



0V> 



showing that it T w interpolates between R§, guaranteeing a good UV behaviour, and R 1 , 
well behaved enough in the infrared. By transposition, using we derive 



T w f{x) 



\a\=k 
\k+l 



fk+d 



'X 

w 1 — 

t 



+ +1) & 

\/3\=k+l 

This is the central formula of this paper. 



t k+d+i ■< \ t 



X 



x 



(6) 



2.2 On the auxiliary function 

It is important to realize what is (and is not) required of the weight function w, apart from a 
good behaviour at the origin: in view of the smoothness and good properties of / away from 
the origin, we have a lot of leeway, and, especially, w does not have to be a test function, 
nor to possess compact support. Basically, what is needed is that w decay at infinity in the 
weak sense that it sport momenta of sufficiently high order. 

We formalize this assertion for greater clarity. First, one says that the distribution / is 
of order \x\ l (with I not a negative integer) at infinity, in the Cesaro sense, if there exists a 
natural number N, and a primitive Jn of / of order N, such that Jn is locally integrable for 
\x\ large and the relation 

f N (x) = 0(\x\ N+l ) 

as \x\ t 00 holds in the ordinary sense. Now, for any real constant 7, the space /C 7 is 
formed by those smooth functions <p such that d a (j)(x) = 0(|x| 7 ~' a ') as \x\ f 00, for each \a\. 
A topology for /C 7 is generated by the obvious family of seminorms, and the space K, is 
defined as the inductive limit of the spaces /C 7 as 7 j 00. Consider now the dual space K,' of 
distributions. The following are equivalent E7| : 
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/e/C'. 
/ satisfies 

f(x) = o(\x\ 

in the Cesaro sense as \x\ | oo. 



— OON 



There exist constants /x Q such that 

/(Ax) ~ £ 



// a 5 (a) (z) 

a>0 



Al»l+ 1 



in the sense of distributions, as A | oo. 



• All the moments (f(x), x a ) exist in the sense of Cesaro summability of integrals (they 
coincide with the aforementioned constants 

Any element of KJ which is regular and takes the value 1 at zero qualifies as a weight 
"function". For instance, one can take for w an exponential function e tqx , with g / 0. 
This vanishes at oo to all orders, in the Cesaro sense, and so it is a perfectly good infrared 
problem-buster auxiliary function. The fact that e iqx £ JC', for q ^ 0, means that, outside 
the origin in momentum space, the Fourier transform of elements £ K, can be computed 
by a standard Cesaro evaluation 

(j)(q) = (exp(iqx),<p(x)). 

Of course, for this auxiliary function the original equation (^) no longer applies, since it has 
no vanishing derivatives at the origin. But (|j) can be replaced by the more general 



W w (f>(x) := <j>(x) - w(x) V ^(^j 



(a) 

(0). (7) 



o<H<fc 



This was seen, at the heuristic level, by Prange p8[ ; see the discussion on the BPHZ formal- 



ism in subsection 5.3, where the "Cesaro philosophy" comes into its own. 

These observations are all the more pertinent because the contrary prejudice is still 
widespread. For instance, the worthy thesis PD| , despite coming on the footsteps of yet 
unfortunately exhibits it; on its page 30: "... [the exponential] function is not allowed in 
the P^-operation because it does not have compact support." Of course it is allowed: then 
the Fourier transformed subtraction W w of Epstein and Glaser becomes exactly the standard 
BPHZ subtraction, around momentum p = q. What T w becomes will be revealed later. 

2.3 Properties of the T-projector 

Consider now the functional variation of the renormalized amplitudes with respect to w. 
One has 
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for T w , and similarly for W w , by definition of functional derivative. It is practical to write 
now 

, i 5 (Q) 
al 

for this combination is going to appear with alarming frequency. From (|4|) we would obtain 

^W„/[.] = -£/[*"•] (8) 

|a|<fe 

whereas (||) yields 

^/[•] = -£/[*".]*», 

\a\=k 

independently of w in both cases. Malgrange's theorem says that different renormalizations 
of a primitively divergent graph differ by terms proportional to the delta function and all its 
derivatives S^ a \ up to |a| = k. Thus there is no canonical way to construct the renormalized 
amplitudes, the inherent ambiguity being represented by the undetermined coefficients of 
the <5's, describing how the chosen extension acts on the finite codimension space of test 
functions not vanishing in a neighborhood of 0. 

There is, however, a more "natural" way — in which the ambiguity is reduced to terms 
in the higher-order derivatives of 5, exclusively. This is guaranteed by our choice of T w . 

In practice, one works with appropriate 1-parameter (or few-parameter) families of aux- 
iliary functions, big enough to be flexible, small enough to be manageable. Recall than in 
QFT, with c = h = 1, the physical dimension of length is inverse mass. Let then the variable 
H have the dimension of mass. We consider the change in T w f when the variable w changes 
from w = w(fj,x) to w((fi + <5/i)x), which introduces the Jacobian jjt = dw ^ x > ; yielding 

|- T w(flx) f = x v d v w{vx)x«) 5 a . (9) 

" \a\=k 

Here we have assumed that / has no previous dependence on fi. 

Enter now the (rotation-invariant) choice w^(x) := H(fi~ 1 — \x\), where H is the Heaviside 
step function: it not only recommends itself for its simplicity, but it turns out to play a central 
theoretical role. The parameter /i corresponds in our context to 't Hooft's energy scale in 
dimensional regularization — see subsection 5.1; the limits n I and \x | oo correspond to 
the case w — 1 and respectively to the "principal value" of /; in general they will not exist. 

Write T M / for the corresponding renormalizations. With the help of @, one obtains 

' |a|=/c 

For / homogeneous (of order — d — k as it happens), the expression is actually independent 
of fi, the coefficients of the 5 a being 



fx a = / f\x\x a , (11) 

\x\=l J\x\=A 
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with \a\ = k and any A > 0. Note that similar extra terms, with \a\ < k, coming out of the 
formulae (|8|) would indeed be //-dependent. 

Compute the T M in the massless (homogeneous) case, whereupon one can pull / out of 
the integral sign. We get 



T,f(x) = R fe { £ d« (1 - (1 - »\x\) k ) 



\a\=k 



+ {k+i) 9(3 

|/3|=fc+l 



X 13 f{x) 



li\x\ 



dt 



t 



(12) 



Formula ( |i~2|) is simpler than it looks: because of our previous remark on (|10|) , all the /> 
polynomial terms in the previous expression for T^f must cancel. Let us then denote, for 
k > 1, 



i=i 



I 



1+1 'k 
I 



1 (k 

2 V2 



+ ■■•-- 



(13) 



At least for < 1, the expression for T„/ becomes 



! (A; + 1) ^ ^ 

|/3|=fc+l 



x 13 f{x) 



(log /x I a; I + 



By performing the derivative with respect to log fi directly on this formula, one obtains 
in the bargain interesting formulae for distribution theory. Namely, for any / homogeneous 
of degree —d — k: 



\p\=k+i 



x 13 f{x) 



E 

\a\=k 



x =1 



fx a )S a (x). 



(14) 



Thus the final result is 



T,f{x) = {-)\k + l) & 

|/3|=fe+l 



X 13 f{x) 



\ogfi\x\ 



+ H k J2c a S a (x) 

\a\=k 



(15) 



with the c a given by (pTTJ) . The resulting expression is actually valid for all x: away from the 
origin it reduces to f(x), as it should. It is in the spirit of differential renormalization, since 
/ is renormalized as the distributional derivative of a regular object that coincides with / 
away from the singularity. 

Let us put equation (13) to work at once. When performing multiplicative renormaliza- 
tion in the causal theory ||, the relevant property of a renormalized amplitude turns out 
to be its dilatation or scaling behaviour. This is not surprising in view of the form of our 
integral equation (|). Now, a further consequence of the choice of operator T w is that it 
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modifies the original homogeneity in a minimal way. Had we stuck to W w (using (|i4|)), the 
relatively complicated form 

WJ(X-) = \- k - d (w^f + log A ° a 6a + E aa (^~ lal ~ !) 5 °) > 
for some a a , would ensue; whereas for T M from fll5|) one obtains 

L/W(Az) := T„/(Az) = \- k - d ([f] R>li {x) +logA E c*r(x)Y 

^ |a|=fc ' 

or 

:= (-* - d)[/WC0 + E 

|a|=fc 

where i? denotes the Euler operator X]|/3|=i x/3 ^- F rom now on , the notations [f]^ and 
7)j/ will be interchangeably used; when the dependence on \i need not be emphasized, we 
can write \J\r instead. 

We invite the reader to prove flTE|) directly, by reworking the argument used for the case 



(16) 



k = in J?]], pp. 307-308]. 

Let us record here the obvious fact that when we employ two different prescriptions 
compatible with our renormalization scheme, the difference in the results is given by 

[/k = [/k + E^> ( 17 ) 

\(3\=k 

for some constants . "Different prescriptions" could mean that we use different weight 
functions, or different choices of renormalization conditions, or simply different values of the 
parameter /i. In this respect, we note that the Fourier transforms of our real space renormal- 
ized amplitudes (see the end of Section 5) do not obey standard renormalization conditions 
on momentum space, so in particular contexts they might be in need of modification for this 
purpose. 



2.4 Supplementary remarks 

Following |32| , with a difference by a factor of 2 in the definition, we conveniently formalize 
the first order operator from the (larger) space 5{(lR d ) of continuous linear functionals on 
Schwartz functions that vanish at the origin, appearing in the foregoing: 

S .= ^xK 

101=1 

It "regularizes" S[, mapping it onto S'. Clearly, on tempered distributions S = d + E. 
Therefore S kills any homogeneous distribution of order —d, like 5 — but not the homogeneous 
Junctionals] From fll5|), it is clear that (the analog in £-space of) S captures the Wodzicki 
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residue density in the theory of pseudodifferential operators — see for instance the discussion 
in Chapter 7 of ]31|. Analogously, define 



s k+1 ■.= (k + iy. J2 ~ir- 

This sends onto S' the space of functionals iS^, +1 (M d ), dual of the space of Schwartz functions 
which vanish at the origin up to order k + 1. On well-behaved distributions, it is equivalent 
to (S + k) ■ ■ ■ (S + 1)S: we think of Sk+i as an ordered power of Si, with the coordinate 
multiplications remaining to the right of the differential operators. For massless models, our 
formulae come close to simply iterating S in this way, as done in P^| . 

For theoretical purposes, it should be kept in mind that the T-operator remains in the 
general framework of the Epstein-Glaser theory: after all, one could always find weight 
functions such that W w ^) is identical to T M . In particular, the singular order of T^f 
is the same as that of / |14 |. However, we see already that a finer classification needs to be 
introduced. 

Definition 1. A distribution / = f\ is called associate homogeneous of order 1 and degree a 
when there exists a homogeneous distribution f Q of degree a such that 

f 1 (Xx) = X a (f 1 (x)+g(X)f (x)), 

for some function g(X). It is readily seen that only the logarithm function can foot the bill 
for g. Then, a distribution f n is called associate homogeneous of order n and degree a when 
there exists an associate homogeneous distribution f n _i of order n — 1 and degree a such 
that 

f n (Xx) = X a (f n (x) + log A /„_!(*)). 

Clearly the renormalization of primitively divergent graphs in massless theories, using T M , 
gives rise to associate homogeneous distributions of order 1. To pass associate homogeneous 
distributions of order 1 thru the same machine (|T2D, in order to obtain a renormalized closed 
expression, is routine. Assuming no prior dependence of fi on /i, one would have to use now 



\a\=k 

instead of ([H]), to simplify the output of (|). We omit the straightforward details. Of course, 
if g = [^r-h for / primitive, then T^g = g; but diagrams with a renormalized subdivergence 
provide less trivial examples. The complete renormalization of massless theories gives rise 
exclusively to associate homogeneous distributions @. 

As remarked in on the finite dimensional vector space spanned by an homogeneous 
distribution and its associates up to a given order, the Euler operator takes the Jordan 
normal form. The Sk+i operators are nilpotent on that space, when a = —d — k. For 
example, 

1 



S 
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and therefore 5' 2 [^] i j = 0. In consonance with this, equation (|T5| ) remains valid for / 
primitively divergent of the associated homogeneous type. 

The main theoretical development of the T-operator closes here. We have found two in- 
stances in the QFT literature of the improved causal method advocated by us (0): the "spe- 
cial R operation" introduced in the outstanding (apparently unpublished) manuscript |3~| ; 
and, unwittingly, "natural" renormalization [KJ — see subsection 5.2. The scaling properties 
are discussed in [|34 |; but no explicit formulae for the renormalized amplitudes are given 



there. We turn to matters of illustration and comparison. 



3 Some examples 

We compute the simplest primitive diagrams relevant for the four-point function of 0| theory 
— for these logarithmically divergent graphs, there is of course no difference between T w and 
W w . The following notation will be used in the sequel: 



n 



d,m 



2m 



x =1 



2r(m + i)vr( d - 1 )/ 2 
T(m + p) 



where i labels any component. In particular, f^o =: fid, the area of the sphere in dimen- 
sion d. The quotients of the Qd,m are rational: 



n 



d,m 



(2m- V 



n d (2m + d-2)(2m + d-A)---d 
The "propagator" , given by the formula 

\x\ 2 ~ d 

when d ^ 2, and by D F (x) = (log |x|)/fi 2 when d = 2, is simply the Green function for the 
Laplace equation: 

ADp(x) = -S(x). 

Consider the "fish" diagram in ip\ theory, giving the first correction to the four-point function. 
The corresponding amplitude is proportional to x~ 4 . On using (|I5|) with k — 0, 



" 1 " 




x 13 \og(/u,\x\) 


.2*. 


R 


X 4 



Si 



log(^\x\) 



X' 1 



Next we look at fij- (T^Dj^j . By direct computation, on the one hand, 



d 



and on the other, according to (|10D : 



" 1 " 


= d fi 


V 


2 A (x 2 ) 


_x 1 _ 


R 







fi 4 5(x) 



d 

"an 



x> 



c°S(x), with residue c =(—,S(/i 1 



R 



X \x\ 



x> 
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which serves as a check. In this case, as /i varies from to oo, all possible renormalizations 
of Dp are obtained. One can as well directly check here the equation: 



x 



(19) 



for f(x) = l/x 4 . Equation (|T9|) contains the single most important information about the 
fish graph and is essential for the treatment of diagrams in which it appears as one of the 
subdivergences [0]. 

Next among the primitive diagrams relevant for the vertex correction comes the tetrahe- 
dron (also called the "open envelope") diagram. In spite of appearances, it is a three-loop 
graph, as one of the circuits depends on the others; it is the lowest-order diagram in (fif 
theory with the full structure of the four-point function. The unrenormalized amplitude 



pnmV 



is of the form 



f(x,y,z) 



x 2 y 2 z 2 (x — y) 2 (y — z) 2 (x — z) 2 ' 

which is logarithmically divergent overall. A funny thing about this diagram is that the 
amplitude for it looks exactly the same in momentum space — see for instance [35 . Denote 
by s the collective variable (x,y,z). Then 



IfU = s t 



log(/i|s| 



x 2 y 2 z 2 [x — y) 2 (y — z) 2 [x — z) 2 



Again, the most important information from the diagram concerns its dilatation proper- 
ties. Proceeding as above, we obtain 



where, for any A > 0, 



^[f] R =:I6(x) 



Is I ds 



2 =A x 2 y 2 z 2 (x — y) 2 (y — z) 2 (x — z)^ 



This integral is computable with moderate effort. First, one rescales variables: y = \x\u and 



\x\v, to obtain 



d 



d 4 u d 4 v 



u 2 v 2 (x/\x\ — u) 2 (u — v ) 2 (x/\x\ — v)' 



The calculation is then carried out by means of ultraspherical polynomial 0, [37| techniques. 
We recall that these polynomials are defined from 



'1 — 2xr + r 



2\-n 



k=0 



x)r 
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for r < 1. There follows an expansion for powers of the propagator 



1 



oo 

|y| fc=0 



(x — ?/) 2n 

if, for instance, |x| < \y\. Using their orthogonality relation 

xu \ ( xv 



d 



S' 



\x\\u\ 



C? 



\x\\v\ 



°kl^k 



UV 



n Q r . 



u\ \v I I k + n 



to perform the angular integrals (in our case n — 1, I — 3), we obtain I as the sum of 
six radial integrals, corresponding to regions like \u\ < \v\ < 1, and so forth. Each one is 
equivalent to 27r 6 C(3). This yields finally the residue 127r 6 £(3) — the geometrical factor fl 4 
is always present. In consequence, now 

TJiXx) = \- 12 [T,f{x) + 12A(3)logA5(x)]. 

This is the first diagram which has a nontrivial topology, from the knot theory viewpoint, 
and thus the appearance of a £- value is expected [pSfl . 

Consider now the two-loop "setting sun" diagram that contributes to the two-point func- 
tion in iff theory; it will prove instructive. One has to renormalize 1/x 6 , and the singular 
order is 2. Off (1151) we read that 



1 

x 6 



J R 



log(ii\x\ 



X v 



37r2 a r/ \ 

+ —A5{x). 



(20) 



" 1 " 


(Ax) = A~ 6 ( 


" 1 " 






_x 6 _ 


R V 


-X 6 - 



Clearly, our formulae come rather close to simply iterating the operator S, as done in [32 
The last term obviously does not make a difference for the dilatation properties; but we shall 
soon strengthen the case for not dropping it. One has 

- ^logA A5(x] 



The reader may check that using W w instead of T w would bring to (^) the extra term 
7T 2 /i 2 <5(x), with an unwelcome /i-power dependence. As we know, this complicates the dilata- 
tion properties for the diagram. The terms polynomially dependent on /z are like the "junk 
DNA" of the Epstein-Glaser formalism, as they carry no useful information on the residues 
ofQFTi. 

More generally, for quadratic divergences (such as also appear in the first (two-vertex) 
contribution to the two-point function of the and y?g theories), one constructs the extension 



2 3 



log/i|x| 



X 



d+2 



(21) 



and 



x 



)=A 



-d-2 



X 



-d-2 



] R + Yd logXA5{x] 
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4 The comparison with the mathematical literature 



4.1 On the real line 

Now we must muster support for the choice of T w and T M . For the basics of distribution 
theory, we recommend J5H]. For concrete computations, a good place to start is the treatment 
by Hormander in Section 3.2 of |39|], of the extension problem for the distributions 

f(x) = x+\ xZ\ M~\ |x| _/ sgn(x), x~ l 

on the real line. Of course, these are not independent: x_ is just the reflection of x + with 
respect to the origin, then \x\~ l = x+ l + xZ l , and so on; note that x~ x is just the ordinary 
Cauchy principal value of 1/x. 

On our side, for instance, xf(x) = H(x) for f(x) = x^ 1 ; xf(x) = sgn(x) for f(x) = 
and so on. Then our formulae flT5|) , for I odd, give 

or, say, 

and, for I even, simply 

W-']* = Sl? 1 ^. (22 > 

Hormander invokes the natural method of analytic continuation of x\, with z complex, 
plus residue subtraction at the simple poles at the negative integers. Our formulae coincide 
with Hormander's — see for example his (3.2.5) — provided that (a) we take fi = 1; and (b) 
.f/fc defined in ( |13D equals (as anticipated in the notation) the sum of the first k terms of the 
harmonic series! This turns out to be the case, although the proof, that the curious reader 
can find in 0, Ch. 6], is not quite straightforward. Thus we understand that in fllBD and 
similar formulae just means X^=i Vj- 

Encouraged by this indication of being on the right track, we take a closer look at the 
analytic continuation method. The point is that the function z i— > J °° x z (j)(x) dx for 9ftz > — 1 
is analytic, its differential being dz J °° x z \ogx <j)(x) dx. Let us now consider the analytic 
continuation definition for x z +1 where for simplicity we first take —2 <dtz < —1. One gleans 

POO 

(x+, 4>) = / x z Rq4>(x) dx. 
Jo 

We recall the proof of this: 

: = <7TT^^ + ^> = -ITT 'ta/°°^ + V'W^. 

A simple integration by parts, taking v = x z+1 and u = <f>(x) — (f)(0), completes the argument. 



[l-l)\dx l 



(sgn(z) log(//| x\) + 2Hi- X 8 



i-i 
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Iterating the procedure, one obtains 

poo 

(x z + ,(p) = / x z R l u l <p(x)dx 
Jo 



for —I — 1 < $lz < — /, with / a positive integer. At z — —I, however, this formula fails 
because of the attendant infrared problem. Let us then compute the first two terms of the 
Laurent development of x z : in view of 

M" 1 r°° . _ /A('- 1 )(0)^-(*+0 



poo 

(x z + ,(f))= / x z R l Q - 1 (f)(x)dx+ / R l Q - 2 (t){x)dx 

J0 Ja- 1 



(Z-l)!(z + Z)' 



the pole part is isolated. Therefore 



lim 

z-v-J 



7^(0 - ^(x) log/i. 



(i-l)!(* + 

Hormander goes on to consider Hadamard's /imie pari: that is, for x+, one studies 



i: 



x z (p(x) dx, 



where is always a test function, for any z G C, and discards the multiples of powers e~ e , 
for nonvanishing with 310 > 0, and the multiples of loge. He proves that this finite part 
coincides with the result of the analytic continuation method. 

We do not need to review his proof, as we can show directly the identity of our results 
with finite part, by the following trick: 

/[i- 1 roo 
x~ l R l ' 1 (f)(x)dx + J x- l R l 2 <p(x)dx. 

Then, as e J, 0, the two last terms give rise to the T^(l/x l ) renormalization and the surviving 

finite terms cancel, except for the expected contribution — ^^yp log/i, coming from the 
first sum. 

Denote the finite part of x+ l by Pf ^j^, where Pf stands for pseudofunction (or for partie 
finite, according to taste). In summary, we have proved: 

Proposition 1. On the real line, the T -operator leads to a one-parameter generalization of 
the finite part and analytic continuation extensions, to wit, 

[x~ l ] R := T,(x- 1 ) = Pf + 5 l '\x) log/i. 

Ob 
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This generalization is in the nature of things. Actually, the finite part and analytic 
continuation methods are not nearly as uniquely defined as some treatments make them 
appear. For instance, at the negative integers the definition of the finite part of x z changes 
if we substitute Ae for e; and, analogously, one can slip in a dimensionful scale in analytical 
prolongation formulae. The added flexibility of the choice of \i is convenient. 

We parenthetically observe that the nonhomogeneity of T^, and then of Pf, is directly 
related to the presence of logarithmic terms in the asymptotic expansion for the heat kernels 
of elliptic pseudodifferential operators f26fl . 

Finally, we remark that the Laurent development for xi| 2= _^ continues: 



^(O) 



+ Pf (xl') + J Pf (x~ l log x+) + E - Pf {x~ l log 2 x+) + • • • (23) 



2 



e(l- 1)! ' -~ v "+ y 1 2 " v ~ + ' 3! 

with e := z + I and the obvious definition for 

Pi(x- l \og m x + )=[x- l \og m x + ] R ^ =1 . 

4.2 Dimensional reduction 

The phrase "dimensional reduction" is used in the sense of ordinary calculus, it does not refer 
here to the method of renormalization of the same name. The reader may have wondered 
why we spend so much time on elementary distributions on R. The reason, as it turns out, is 
that an understanding of the 1-dimensional case is all that is needed for the renormalization 
of |x| _,i_A: , for any k and in any dimension d; thus covering the basic needs of Euclidean field 
theory. For instance, one can define [^ _4 ] i? on R 4 from knowledge of X+ 1 on R. 

Denote r := \x\ and let f(r) be an amplitude on R d , depending only on the radial 
coordinate, in need of renormalization. We are ready now to simplify flT5| ) by a method that 
generalizes Proposition 1 to any number of dimensions. 

Given an arbitrary test function <ft, consider its projection onto the radial-sum-values 
function <p i— ► Pcj) given by 



P<p(r) := / (j)(ry). 
J\v\=i 

We compute the derivatives of Pcj) at the origin: 

(P0)(2m+i)( O ) = o and 
(P0)( 2 ™) (O ) = ^ m A m 0(O). 
To prove this, whenever all the /3's, and thus n, are even, use 

2 rf ) Y( l3n+1 
v?---v p n = ' ' ' ' x " 

x =o J\y\=i 



\/3\=n P - x=0J\y\=l 



r(^) 

in consonance with (]T8|); the integral vanishes otherwise. Note that Pct> can be considered as 
an even function defined on the whole real line. Then, whenever the integrals make sense, 

(fi^^ix^^ifiry-^P^r))^, 
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which in particular means that extension rules for H(r)f(r) on K give extension rules for 
f[r)r d ~ l on W 1 . This we call dimensional reduction. 

Before proceeding, let us put the examined real line extensions in perspective, by inves- 
tigating how satisfactory our results are from a general standpoint, and whether alternative 
renormalizations with better properties might exist. Note first, from fllBl) : 

pf EM. = 4-(H(x) log Id). 
x dx 

For z not a negative integer, the property 

xx z + = x z + l (24) 
obtains; and excluding z = as well, we have 

4-xi = zx z ~ 1 - (25) 
dx 

One can examine how the negative integer power candidates fare in respect of these two 
criteria: of course, except for x~ l , which keeps all the good properties, homogeneity is 
irretrievably lost. 

Actually, it is x+ l that we need. One could define a renormalization diff of x+ l simply 

by 

so [^+ 1 ] diff = Pf ; and automatically the second ( p5f) of the requirements 

fa [ x + ] diff = ~l [ x + ] diff 

would be fulfilled. This would be "differential renormalization" in a nutshell. It differs from 
the other extensions studied so far: from our previous results, 

Ma = [x+] R + (-)'(^-i + log/x) S l -\x). 
On the other hand, it is seen that 



so that Tp does not fulfill that second requirement; but in exchange, it does fulfill the first 
one (0): 

There is no extension of x°l for which both requirements simultaneously hold. 

It looks as if we are faced with a choice between [-]diflf and T M (-) — which is essentially 
Pf(-) — each one with its attractive feature. But the situation is in truth not symmetrical: 
in higher dimensional spaces the analogue of the first requirement can be generalized to the 
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renormalization of \x\ ; whereas the analog of the second then cannot be made to work 
— have a sneak preview at (|26|). 



Estrada and Kanwal define then, for k > [}4T 

1 



Pf 



~d+k 



1 



„d+k 



diff 



0(x) 



Pf 



diff 



,P0(r 
P0(r) 



In view of (^2j), the case odd is very easy, and then all the definitions coincide: 



Pf 



~d+k 



T, 



V \ r d+k 



„d+k 



diff 



-d-fc' 



the function r z having a removable singularity at — d — k. However, in most instances in 
QFT k happens to be even, so we concentrate on this case. We are not in need of new 
definitions. By going through the motions of changing to radial plus polar coordinates and 
back, one checks that, assuming a spherically symmetric weight function w, the evaluation 
(T w f(r), 4>{x)) is equal to 



/(r), -0(0) 



A0(O) 
2\d 



r 2 — ■ ■ ■ — w(r) 



O d , m A m 0(O) 



,2m 



the right hand side being invariant under T w . This was perhaps clear from the beginning, 
from symmetry considerations. It means in particular that the different putative definitions 
of on M d obtained from T M on the real line all coincide with the original definition, that 
is: 

Proposition 2. The operators are mutually consistent under dimensional reduction. 
Moreover, 



r^T^r 



-d—2m\ 



rp ^-d-2m+2q^ 



follows, by using the easy identity 



r 2 A m 5(x) = 2m(2m + d- 2)£ m ~ 1 5(x). 



Therefore, it is now clear that 



-d—2m\ 



Pf(r 



-d—2m\ 



Q dtm A m 6(x) 
^ (2m)! 



log/i. 



It remains to compute the derivatives. A powerful technique, based on "truncated regu- 
larization" and calculation of the derivatives across surface jumps, was developed and clearly 
explained in pfl|| . It is rather obvious that for k — d odd the "naive" derivation formulae 
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(see right below) will apply. Whereas for k — d = 2m even, they obtain extra delta function 
terms; in particular for the powers of the Laplacian 



1 



d+2m 



(d + 2m + 2n - 2) ■ ■ • (d + 2m + 2)(d + 2m)(2m + 2) 



J diff 



x (2m + 2n) 



1 



r 



d+2m+2n 



+ 



n 



d,m 



diff 



A n S(x) 



(2m)! 2m + 21 - 1' 



(26) 



The first term is what we termed the "naive" formula. 

Estrada and Kanwal do not explicitly give the powers of A for finite part. But from (|26|) 
is a simple task to compute 



A n 



„d+2m 



(d + 2m + 2n - 2) ■ ■ ■ (d + 2m + 2)(d + 2m)(2m + 2) 



J R 



x (2m + 2n) 



r 



d+2m+2n 



d,m 



(4(m + Q + d - 2) A n S(x) 

(2m)! 2(m + I) (2m + 21 + d - 2) 



No one seems to have computed explicitly the distributional derivatives of the 
Pf(af^ log m a; + ) and the correspondingly defined Pf (r~ l log" 1 r), although they might be quite 
helpful for Euclidean QFT on configuration space. 

We next enterprise to tackle a comparison with methods of renormalization in real space 
in the physical literature. Of those there are not many: it needs to be said that the flame- 



keepers of the Epstein-Glaser method [43[] actually work in momentum space (using dis- 
persion relation techniques). Euclidean configuration space dimensional regularization, on 



the other hand, starting from |37j] , evolved into a powerful calculational tool in the eighties. 
With the advent of "differential renormalization" in the nineties, regularization-free co- 
ordinate space techniques came into their own: they are the natural "market competitors" 
for the ideas presented here. 

We deal first with dimensional regularization. 



5 Comparison with the QFT literature 

5.1 Dimensional regularization and "minimal subtraction" 

Dimensional regularization on real space, for primitively divergent diagrams, can be identified 
with analytic continuation. To get the basic idea, it is perhaps convenient to perform first a 
couple of blind calculations. Start from the identity 

A\- d+e ^^S x (\x\- d+e ). 
e 

Then, expanding in e, on use of (|ToD, it follows that 

e \xr 
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The first term is a typical infinite (as 1/e) counterterm of the dimensionally regularized 
theory. The order of the delta function derivative, in this case, tells us that we are dealing 
with a logarithmic divergence. The coefficient Q n of the counterterm, or QFT residue, 
coincides with our scaling coefficient of Section 2. The second term is precisely our 
renormalized expression. 

Let us go to quadratic divergences. A brute-force computation establishes for them the 
differential identity 

^\~ d - 2+e = ^ £ , i 2 , S 3 (\x\ M ). (27) 



On the other hand, from (IT 



2e(l-fe+ie 2 



S 3 {\x\- d - 2 )=^A5(x). (28) 



d 

Performing in (|2T|) the expansion with respect to e, this yields 

-d-2+e _ Ar/ \ | 1 c flog(l\x\\ 3£l d 



That is, 



A pattern has emerged: as before, there is a unique counterterm in -; the residue coincides 
with our scaling coefficient; the order of the delta function derivative reminds us of the order 
of the divergence we are dealing with; and the "constant" regular term is precisely [l/|x| d+2 ].R 
constructed in (pip according to our renormalization scheme. 

The correspondence between the two schemes, at the present level, is absolute and 
straightforward. It is then a foregone conclusion that we shall have //-independent residues, 
always coincident with the scaling factors, for the simple poles of l/\x\ d+2m , and that the first 
finite term shall coincide with T M , provided we identify our scale with 't Hooft's universal 
one. This is an immediate consequence of the Laurent development (p3|) , transported to M, d 
by dimensional reduction. In symbols 

S\x\ Mm = ^l^^) + [\*\- d - 2 ] R + 0(e). (29) 

This substantiates the claim that T M effects a kind of minimal subtraction. Let us point out, 
in the same vein, that already in |2£| the analytic continuation of Riemann's zeta function 
was evaluated as the outcome of a quantum field theory-flavoured renormalization process. 

A word of warning is perhaps in order here. Performing the Fourier transform of these 
identities, we do not quite obtain the usual formulae for dimensional regularization in mo- 
mentum space. The nonresemblance is superficial, though, and related to choices of "re- 
normalization prescriptions". The beautiful correspondence is "spoiled" (modified) as well 
for diagrams with sub divergences, because in dimensional regularization contributions will 
come to 0(e°) from the higher terms of the e-expansion, when multiplied by the unavoidable 
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singular factors; but, again, the difference is not deep: we show in III how one organizes 
the Laurent expansions with respect to d so as to make the correspondence with the T- 
subtraction transparent. 



Much was made in [53], and rightly so, of the importance of the perturbative residues 
in the dimensional regularization scheme. Residues for primitive diagrams are the single 
most informative item in QFT. The coefficients of higher order poles are determined by the 



residues — consult the discussion in Now, the appeal of working exclusively with well- 
defined quantities, as we do, would be much diminished if that information were to disappear 
in our approach. But we know it is not lost: it is stored in the scaling properties. 

5.2 Differential renormalization and "natural renormalization" in QFT 

Differential renormalization, in its original form, turns around the following extension of 
l/x A (in M 4 ): 



[l/.V^-iA^. (30) 
4 x A 



At present, two main schools of differential renormalization seem still in vogue: the original 
and more popular "(constrained) differential renormalization" of the Spanish school — see 
for instance [fJBJ — and the "Russian school" — inaugurated in |52| . This second method, as 
already reported, reduces to systematic use of the operators Sk+i, i.e., to our formulae (p~5|) 
without the delta terms. Whereas the first school has its forte in concrete 1-loop calculations 
for realistic theories, assuming compatibility of differentiation with renormalization, the 
second initially stressed the development of global renormalization formulae for diagrams 
with subdivergences, and the compatibility of Bogoliubov's rules with renormalization. 

Hereafter, we refer mainly to the original version. It proceeded from its mentioned start- 
ing point to the computation of more complicated diagrams by reductions to two-vertex 
diagrams. This involves a bewildering series of tricks, witness more of the ingenuity of the 
inventors than of the soundness of the method. V. gr., the tetrahedron diagram (considered 
already) is rather inelegantly renormalized by the substitution 1/x 2 i— > x 2 [1/x 4 ]r,,fjl- They 
get away with it, in that particular case, because their expression is still not infrared diver- 
gent. But in nonprimitive diagrams infrared infinities may arise in relation with the need 
to integrate the product of propagators over the coordinates of the internal vertices in the 
diagram, and, in general, under the procedures of differential renormalization it is impossible 
to avoid incurring infrared problems M7\ . 



Even for primitively divergent diagrams, differential renormalization is not free of trouble. 



In his extremely interesting paper, Schnetz p5| delivers a critique of differential renormali- 
zation. In elementary fashion, notice that 



xV\og(n 2 x 2 ) 1 



x A 2 

and so 



— 0" 



\og(n 2 x 2 ^ 



X 2 
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This is to say: 



[V^ R -[V*\ FJL = ^)- (31) 

We contend that "our" [l/x 4 ]ij and not [1/x 4 ]r,fjl is the right definition. Of course, one is 
in principle free to add certain delta terms to each individual renormalization and proclaim 
that to be the "right" definition. However, 1/x 4 on M 4 is dimensionally reduced to x+ l on 
M + and because, as already pointed out, differential renormalization of this distribution is 
consistent with [x^ 1 ]^ for fj, — 1, the [-]r,fjl definition is inconsistent with any of the natural 
alternatives we established in the previous subsection. (It would clearly induce back an extra 
S term in the definition of [x+^r on the real line, fully unwelcome in the context.) 

In other words, if we want to make use both of sensible rules of renormalization for the 
radial integral (namely, including differential renormalization at this level) and of Freedman, 
Johnson and Latorre's formulae, we have to relinquish the standard rules of calculus. This 
Schnetz noticed. 

Schnetz proposes instead a "natural renormalization" procedure on M 4 , boiling down to 
the rule 



x 2 



A' 



in+l 



n\(n + l)\ 

x 2n+4 



R 



2 



[8n 2 H n + ^-)A n 6(x), (32) 
n 



whose first instance is precisely the previous equation (plf ). This he found by heuristically 
defining "natural renormalization" as the one that relates renormalization scales at different 
dimensions without changing the definition of ordinary integrals or generating r-dependence 
in the renormalization of r-independent integrals; and by elaborate computations to get rid 
of the angular integrals. 

His calculation is any rate correct, and the results can be read off (for d — 4, m — 0) 
our fl2HD , taking into account (|30| ) and (|31|). We have proved that our operator T w in the 
context just amounts to "natural renormalization". 

Shortly after the inception of the differential renormalization, it naturally occurred to 
some people that a definite relation should exist between it and dimensional regularization. 
However, because of the shortcomings of the former, they landed on formulae both messy 



and incorrect [48]. The reader is invited to compare them with our p 



The more refined version of differential renormalization in 15211, coincides with our for- 



mulae for logarithmic divergences and eludes the main thrust of Schnetz's critique; however, 
we have seen that in general it does not yield the Laurent development of the dimensionally 
regularized theory either. On the other hand, it must be said that the emphasis in |24|, [25| in 
bringing in the Laplacian instead of the less intuitive albeit more fundamental Sk operators 
has welcome aspects, not only because of the enhanced feeling of understanding, but also 
in that it makes the transition to momentum space a trivial affair, as soon as the Fourier 
transform of the (evidently tempered) distribution x~ 2 log(/i 2 x 2 ) is known. 

The trinity of basic definitions in differential renormalization is then replaced by the 
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identities 



r 1 
1 

log fi\x\ 
X 4 



1 log/i|x| 



■1a« 

16 



log 

X 2 



7r 2 <5(x); 

+ |vr 2 A5(x) 



R 



1 log 2 yu|x| + log yt!\x\ 
"I 



+ y^) 



the <5's being absent in standard differential renormalization. In the next Section 6 we shall 
see another demonstration of their importance. 

The kinship of the EG method with differential renormalization a la Smirnov and Zavialov 
was recognized by Prange [28 ; he was stumped for nonlogarithmic divergences, though. 
See [E3| in the same vein. 



5.3 The connection with BPHZ renormalization 

We still have left some chips to cash. We elaborate next the statement that BPHZ subtraction 
has no independent status from Epstein-Glaser, and that the validity of that renormalization 
method is just a corollary of the latter. This involves just a two-line proof. 

The Fourier transforms of the causally renormalized amplitudes exist at least in the 
sense of tempered distributions. They are in fact rather regular. Taking Fourier transforms 
is tantamount to replacing the test function by an exponential, which, according to the 
Cesaro theory of [^6|, [27]] , can preclude smoothness of the momentum space amplitude only 
at the origin. The appearance of an (integrable) singularity at p = is physically expected 
in a theory of massless particles. 

Let us fix our conventions. We define Fourier and inverse Fourier transforms by 

F[0](p) :=#(p):= J ^e^(x), 

and 

F-^Kp) :=kp) ■= J J^^m, 
respectively. It follows that 

(x^yip) = H) M 0"to, 

where \i denotes a multiindex; so that, in particular, 

(x»y(p) = (-iy(2n) d / 2 5^\p). 

Also, 

d"<j)(0) = (-zf(27r)- d / 2 (^,0). (33) 
From this and the following consequence of (Q): 

(F[/],F- 1 [0]) = (F[/],F- 1 [ J R^]), 
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there follows at once 

F[f](p) = R k F[f](p). 

This is nothing but the BPHZ subtraction rule in momentum space. 
We hasten to add: 

• An expression such as F[f] is not a priori meaningless: it is a well defined functional 
on the linear subspace of Schwartz functions whose first momenta f p a (p(p) d d p up 
to order k + 1 happen to vanish. (This is the counterpart FSk+i, according to ([33]), of 
the distributions on real space acting on test functions vanishing up to order k + 1 at 
the origin.) 

• Moreover, explicit expressions for these functionals on the external variables are given 
precisely by the unrenormalized momentum space amplitudes! 

This circumstance constitutes the (deceptive) advantage of the BPHZ formalism for 
renormalization. We say deceptive because — as persuasively argued in |53J — the 



BPHZ method makes no effective use of the recursive properties of renormalization 
(paper III) and then, when using it, prodigious amounts of energy must go into proving 
convergence of, and/or computing, the (rather horrendous) resulting integrals, into 
showing that the Minkowskian counterparts define bona fide distributions . . . Much 
more natural to remain on the nutritious ground of distribution theory on real space, 
throughout. But this has never been done. 

• Also, the <9 M -F[/](0) for < k exist for massive theories. 

For zero-mass models, the basic BPHZ scheme runs into trouble; this is due naturally 
to the failure of 9^/(0) to exist for = k, on account of the infrared problem. Now, one 
can perform subtraction at some external momentum q ^ 0, providing a mass scale. This 
is just the Fourier-mirrored version of standard EG renormalization, with weight function 
e~ iqx ; one only has to remember to use (0) instead of @. 

It is patent, though, that this last subtraction is quite awkward in practice, and will 
introduce in the Minkowskian context a noncovariance which must be compensated by further 
subtractions. This prompted Lowenstein and Zimmermann to introduce their "soft mass 
insertions" Which amounts to an epicycle too many. 



In the light of the approach advocated in this paper, there exist several simpler and more 
physical strategies. 

• One strategy is to recruit our basic formula @ in momentum space 

d^F[fw\{p) p^ 



u! 

\n\=k p 



Fim-jt'Fim-Yl 

Still with w(x) = e~ igx , this leads at once to 
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Note that the difference between two of these recipes is polynomial in p M , with \fi\ = k 
only, as it should. This can be more easily corrected for Lorentz covariance, should the 
need arise 0. 

A second method is to exploit homogeneity in adapting our recipes for direct use in 
momentum space, in the spirit of [|51| and [j5"2X 



A third one is to perform Fourier analysis on our previous results. One has 

d 4 x _ IPX logQuM) J_ (C\p\\ 
(2vr) 2 x 2 p 2 g \2/j )' 

where C := e 7 ~ 1.781072 . . . with 7 the Euler-Mascheroni constant. Then, from 
for instance for the "fish" diagram in the ip\ model: 

[i7?k(p) = |[i-iog(cv/v)], 

and more generally: 



_____ 2k 

[l/x 2 ^] R (p) ' ' ' 



2 log M - V(k + 1) - V(k + 2) 
2/i 



A k+1 k\(k + l)\ 

where := d/dx(logT(x)) has been invoked, and we recall that 

*(n) = - 1 + H n _ l . 
For the setting sun diagram in the tp\ model, in particular: 

m^t-W^ 1 -^- (34) 



16 V 2/i 4 
6 Some examples in massive theories 

The aim of this short section is to dispel any idea that the usefulness of EG-type renormali- 
zation, and in particular of the T-subtraction, is restricted to massless models. The overall 
conclusion, though, is that the massless theory keeps a normative character. Our purposes 
being merely illustrative, we liberally borrow from Schnetz |?5] , Prange [ 2£| and Haagensen 
and Latorre 



The first example is nothing short of spectacular. Suppose we add to our original La- 
grangian for ip\ a mass term |m 2 <^ 2 and treat it as a perturbation, for the calculation of the 
new propagator. Then we would have for Dp(x): 

i_ [ dx > ^ + f dx > dx » ^ 

x 2 J (x — x') 2 x' 2 J (x — x') 2 (x' - x") 2 x" 2 

This "nonrenormalizable" interaction is tractable with our method. We work in momentum 
space, so we just have to consider the renormalization of l/p 2k for k > 1. This is read 
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directly off (|34|), by inverting the roles of p and x, with the proviso that fi gets replaced by 
in order to keep the correct dimensions. Then the result is 



D F (x) 



1 m 



E 



m 2n x 2n 



2 ^4 n ra!(n + l)! 

n=0 V ' 



log 



\l\x\ 



- ^(n + 1) - ^(n + 2] 



On naturally identifying the scale fi = m, one obtains on the nose the exact expansion of 
the exact result 



D F (x) 



m 



\x\ 



Ki(m\x\). 



Here K\ is the modified Bessel function of order 1. Had we kept the original EG subtraction 
with a H(fi — \p\) weight, we would earn a surfeit of terms with extra powers of /i, landing 
in a serious mess. 

It is also interesting to see how well or badly fare the other "competitors" . Differential 
renormalization gives a expression of similar type but with different coefficients: 

m 2n x 2n 



1 m 2 yr-^ m 2n x 2n ( p\x\ \ 
x* + ~2 4-( 4" n\(n + 1)! V ° g + 7' 



2 ^4n n !( n+ l)l 

?1=0 V ' 

To obtain the correct result, it is necessary to substitute a different mass scale //„ for each 
integral and to adjust ad hoc an infinity of such parameters. Dimensional regularization 
(plus "minimal" subtraction of a pole term for each summand but the first) fares slightly 
better, as it "only" misses the \l/(n + 2) terms The distribution-theoretical rationale for 
the success of the "illegal" expansion performed is explained in |[25|| . 

Let us now look at the fish diagram in the massive theory. It is possible to use (fj) instead 
of (P). Make the change of variables: 



\x\ 

s 



dt 



\x\ ds 



One gets, for the renormalized amplitude, 



-S 



m 2 f°° 

— / dssKiOmsY 

% J\x\ 



Now, 



ds sKi{ms) 2 = ■ 
can be easily checked from 

K(s)=K 1 (s) 
The final result is then 



#o(s)#i(s) 



K[(s) = -K (s) - 



-Ki(s] 



KAs) 



A 



2 

— [Kl{m\x\) - Kl{m\x\) + 



K (m 


x\)K 


i(m 


x\) 


rn 


X 
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Had we used (|6]), the upper limit of the integral would become and the result would be 
modified by 



2 

m 



■{Kl{m/n) - K 2 {m/n)Ko{m/n)). 



4:fX 2 X 2 

At the "high energy" limit, as fi j oo and \x\ j 0, this interpolates between the previous 
result and the renormalization in the massless case. 

However, this method becomes cumbersome already for renormalizing Dp. It is conve- 
nient to modify the strategy, and to use in this context differential renormalization, corrected 
in such a way that the known renormalized mass zero limit is kept. This idea succeeds be- 
cause of the good properties of our subtraction with respect to the mass expansion. For 



instance, away from zero |53 





(m 






X 


) 



3 2 

= — (A- 9m z )(A-m 2 )(K (m\x\)Kf(m\x\) +Kn(m\x\)). 
16 

Note the three-particle "threshold" . To this Haagensen and Latorre add a term of the form 

7r 2 2a 

— log A5(x), 

4 7m 

to which, for reasons sufficiently explained, we should add a term of the form ^AS(x). A 
term proportional to 5 (thus a mass correction) is also present. As they indicate, it is better 
fixed by a renormalization prescription. 



7 Conclusion 

We have delivered the missing link of the EG subtraction method to the standard litera- 
ture on extension of distributions. The improved subtraction method sits at the crossroads 
in regard to dimensional regularization in configuration space; differential renormalization; 
"natural" renormalization; and BPHZ renormalization. The discussions in the previous 
sections go a long way to justify the conjecture (made by Connes, and independently by 
Estrada) that Hadamard's finite part theory is in principle enough to deal with quantum 
field theory divergences. To accomplish that feat, however, it must go under the guise of 
the T-projector; this gives the necessary flexibility to deal with complicated diagrams with 
subdivergences [f|. 
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